Search for a light Higgs-hke boson in J/^ radiative decays 
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Using a data sample containing 1.06 x 10* ^p' events collected with the BESIII detector at 
the BEPCII electron-positron collider, we search for a light Higgs-like boson A^ in the process 
tp' — ^ -k'^7v~ J/ij), J/ijj ^ 'yA'^, fi^ fi~ . Particles of this type are predicted in the next- Minimal 

supersymmetric extension of the Standard Model. The HyperCP experiment observed three anoma- 
lous E+ — > events, with fi'^ fi" masses clustered around 214.3±0.5 MeV/c^ . In this analysis, 
we find no evidence for any mass peak between the mass threshold and 3.0 GeV/c^. We set 
90% confidence level upper limits on the product branching fractions for J/ip — >■ jA^, A° — > fi^ fJ-~ 
that range from 4 x 10"'' to 2.1 x 10"^ depending on the mass of A° , for M{A°) < 3.0 GeV/c^. 
Only one event is seen in the mass region below 255 MeY/c^ and this has a fJ,'^ /j.^ mass of 213.3 
MeV/c^ and the product branching fraction upper limit 5 x 10~'. 



PACS numbers: 13.20.Gd, 14.40.Pq, 14.80.Da 
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A light Higgs-like scalar boson is predicted in the 
next-Minimal supersymmetric extension of the Standard 
Model [iHll- Sensitive searches for the are impor- 
tant and can be helpful for a better understanding of the 
Higgs mechanism and supersymmetry. The HyperCP ex- 
periment [3| observed three anomalous — > p^'^ 
events with fJ-^fJ-^ invariant mass clustered around 
214.3 MeV/c^, and these have been interpreted as possi- 
ble evidence for a scalar boson with mass 214.3 MeV/c^ 
that decays into a /x+m" pair @1 . The BaBar 0, , 
CLEO and Belle [l3l experiments have searched for 
A*^ production in various channels and found no evidence 
for a signal. 

For an A° boson with mass below the r-pair threshold 
the decay A*^ — >■ ^jL~ is expected to be more prominent 
than A" — >■ e+e^ because the mass of muon is much 
larger than that of the electron. We use the process ip' 
7r+7r~ J/V', J/tp ^ "yA°, A° — >■ to search for an A° 



111 at the BEPCII electron- 



with the BESIII detector 
positron collider [l2| . 

BEPCII is a double-ring e+e^ collider with a design 
peak luminosity of lO'^'^ cm^^s^^. The BESIII detector 
is based on a large 1 Tesla solenoid magnet and covers 
93% of the total Att solid angle surrounding the e+e" 
collision point with four major detection systems: (1) A 
small-cell, helium-based main drift chamber (MDC) with 
43 layers that provide an average single-hit resolution of 
135 ^m, charged-particle momentum resolution of 0.5% 
at 1 GeV/c, and a dE/dx resolution that is better than 
6%. (2) An electromagnetic calorimeter (EMC) consist- 
ing of 6240 CsI(Tl) crystals configured in a cylindrical 
structure (barrel) and two endcaps. The energy resolu- 
tion for 1.0 GeV 7-rays is 2.5% in the barrel and 5% in 
the endcaps, and the position resolution is 6 mm in the 
barrel and 9 mm in the endcaps. (3) A time-of-flight 
system (TOF) constructed of 5-cm-thick plastic scintil- 
lators, with 176 pieces of 2.4 m long counters arranged 
in a two layer barrel and 96 fan-shaped counters in the 
endcap regions. The barrel (endcap) time resolution of 
80 ps (110 ps) provides 2a K/tt separation for momenta 
up to ~ 1.0 GeV/c. (4)Muon identification is provided 
by 1000 m^ of Resistive Plate Chambers (RPCs) that 
are interspersed in the magnet's iron flirx return (MUC). 
Nine barrel and eight endcap layers provides 2 cm posi- 
tion resolution for penetrating particles. 

The analysis is based on 1.06 x 10^ events collected at 
the peak of the "0' resonance. The number of -0' events 
was determined by counting inclusive hadronic events 
as described in Ref. [l3| with an estimated uncertainty 
of 4%. Monte Carlo (MC) events are simulated with 
the GEANT4 program [1J| and experimentally determined 
resolutions of the wires and counters in the detector. 

For the event selection, we first require two positive 
and two negative charged tracks and at least one good 
photon. We also use ^ identification information, veto 
tt'^s and apply kinematic constraints. 



Candidate photons are energy clusters in the EMC 
that: (1) are within the fiducial region of the EMC 
(I cos 6*..^ I < 0.8 for the barrel and 0.86 < | cos 6*^1 < 0.92 
for the endcaps); (2) are more than 20 degrees away from 
the extrapolated position of the closest charged track; (3) 
have a pulse time that is consistent with being produced 
together with the charged track candidates. 

Charged track candidates are required to originate 

from the interaction point, Vxy =^Jv^~hV^ < 1 cm, 

\Vz\ < 10 cm, where V^, Vy, and are the a;, y, and z 
coordinates of the point of closest approach to the inter- 
action point. The tracks are also required to be within 
the polar angle region | cos0| < 0.93. 

Candidate muons are charged tracks in the active area 
of the barrel MUC (|cos0| < 0.75) with: momentum 
higher than 0.7 GeV/c; energy deposition in the EMC 
between 0.15 GcV and 0.26 GeV; E/p (EMC energy over 
MDC momentum) less than 0.5; at least three associated 
hit layers in the MUC. For tracks in the momentum range 
0.8 GcV/c < p < 1.15 GcV/c, the MUC penetration 
depth is required to be greater than (70p — 40) cm {p in 
GeV/c); for tracks with p > 1.15 GcV/c. the penetration 
depth is required to be more than 41 cm. Tracks with 
momentmir bigger than 0.8 GcV/c are removed if the fit 
to the MUC hits either fails or gives a poor fit result. 

If there are multiple photons with energy larger than 
25 MeV, we reject the event if any pair of these photons 
has an invariant mass within 40 MeV/c^ of to^o. For the 
multi-photon events that remain, the 7 with the highest 
energy is selected as the photon used in the analysis. 
The pair of oppositely charged tracks with recoil mass 
closest to the J/tp mass is assigned as the tt'^ and tt^ 
and the other two tracks as the and /i^. At least 
one of the tracks assigned as a /i is required to satisfy 
the muon candidate criteria. We select events with a 
TT+TT^ recoil mass in the range between 3.092 GeV/c^ 
and 3.102 GeV/c^ and perform a 4-constraint energy- 
momentum conserving kinematic fit using the selected 
7 and four charged tracks. We require < 40 and 
M{n+fj,-) < 3.02 GeV/c2. 

The A'^ particle is expected to have a narrow width. 
Simulations where the width of is set to zero and 
the mass is set at 71 different values that range from 
0.212 GeV/c^ to 3.0 GeV/c^ indicate that the selection 
efficiency varies between 28% and 18%, depending on the 
mass of 

The /i^/i^ mass distribution of selected data events is 
shown in Fig. [BJa). Over the entire mass range, from 
threshold to 3.0 GeV/c^, there is no evident narrow 
peak. Below 255 MeV/c^ there is only one event, with 
fi^H^ invariant mass of 213.3 MeV/c^. The expected 
mass resolution of A° is about 0.2 MeV/c^ for M{A°) at 
213.3 MeV/c^ and the major backgrounds in this region 
come from tp' — > 7r+7r^J/-0, J/ip •yTT^ir^. The ex- 
pected number of background events in the mass region 
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near 213.3 MeV/c^ is about 0.2/MeV/c2; the observa- 
tion of one event in this region is consistent with that 
background leveL 

To set upper Umits on the production rates for different 
masses, we do unbinned maximum hkeUhood fits to the 
fJi^ jjT invariant mass spectrum where the mass of the 
peak is restricted to be within a series of 5 MeV/c^ 
wide bins. In each fit, we use a MC-determined shape 
for the A" signal and for the background shape we use 
a polynomial. Bayesian upper limits on the signal yield 
in each 5 McV/c^ bin are determined. Figure [2] shows a 
typical fit (at 2.43 GeV/c^) and likelihood distribution 
as a function of the signal yield with the 90% C.L. upper 
limit value indicated by an arrow. 

We use different fit ranges, polynomial background 
shapes of different orders, different MC signal shapes for 
different mass values as signal shapes to estimate the 
fit-related systematic error on the signal yield in each 
mass interval. The resulting 90% C.L. upper limits on 
the number of signal events for all the mass intervals are 
shown in Fig. [ijb). 

The systematic errors in the J 
product branching fraction measurement are summarized 
in TablelH these include contributions from tracking, par- 
ticle identification, photon selection, kinematic fit, 7r"*'7r^ 
recoil mass requirement, and tt^ veto. The uncertainty 
of the number of ip' events is 4% O and that of the 

-0' ^ 71 



7^0, AO ^ 11+ II- 



" J /ip branching fraction is 1.2% [15 



TABLE I: The individual contributions to the total relative 
systematic error (%) in the product branching fraction mea- 
surement. 



Source 


Error 


Tracking efficiency 


4.0 


Particle identification 


2.5 


Kinematic fit 


3.2 


7 efliciency 


1.0 


TT+TT^ recoil mass 


1.2 


tt" veto 


2.0 


Number of ip's 


4.0 




1.2 


Total 


7.5 



The uncertainty due to data-MC difference in the 
charged tracking efficiency is 1% per track and added 
linearly. This is determined from samples of J/ip — > pir 
and J/ip — > ppTT+TT^ events. In this analysis, there are 
four charged tracks and the relative systematic error is 
4%. 

The uncertainties due to /i identification (FID) is de- 
termined from studies of a sample of radiative muon pair 
events that contain one photon. The FID efficiency for 
one track is about 65% (during the data taking, the MUC 
end cap detectors were not functioning). The difference 
in fi FID efficiency between data and MC is less than 4% 



per track. Since we require only one muon to satisfy the 
identification criteria, the FID related systematic error is 
less than 2.5%. 

The systematic uncertainty associated with the kine- 
matic fit is determined by applying a similar kinematic 
fit to MC and data samples of tp' tt+tt^ J/tp, J/ifj 
TT+TT^Tr", tt" — !• 77 events. In the event selection for 
this study, if there are more than two candidate 7s we 
use the most energetic 7 together with the one that has 
the best IC fit to 7r° 77. From data-MC differences 
for these events, the systematic error associated with the 
kinematic fit is determined to be 3.2%. 

The systematic error associated with the tt" veto 
is studied with samples of ip' — > jXcj — > j4'4' 
j2{K+K-) and ?/'' ^ tt+tt-J/V', J/V' ^ 7/2(1270), 
72(1270) 7T+7T~ events from both MC simulation and 
data. From ip' IXcJ ~^ I'P't'i we determine the effect 
of the 7r° veto cut on the efficiency. For the second sam- 
ple, we fit the 7r+7r~ mass spectrum to get the number of 
72(1270) — > 7r+7r~ events with and without the 7r° veto 
cut. Data and MC efficiency differences for the tt" veto 
are found to be less than 1.7% from the first channel and 
less than 2.0% from the second channel. We use 2% as 
the systematic error due to the M(77) requirement. 

The systematic error caused by the tt+tt^ recoil mass 
requirement is analyzed with the sample of ip' — > 
TT+'K~ J/ip, J/tp ^ fJ-^ fJ-^ in data and MC simulation. 
From the numbers of events with and without the recoil 
mass requirement we determine data and MC efficiency 
difference to be less than 1.2%. 

The systematic errors are summarized in Table HI As- 
suming the errors from all sources are independent, the 
total error is determined from the quadrature sum to be 
7.5 percent. 

We determine the upper limit on the branching frac- 
tions oi J/4' 7^": from the relation 



B < 



Nsig(UL)/e 



N{ij/) X B{ip' 



T+TT-J/i)) X (1 -cr)' 



(1) 



where Nsig(UL) is the upper limit on the number of sig- 
nal events in each M{fi+ii~) bin after consideration of 
the mass fitting systematic errors shown in Fig. [Ub); e 
is the A°-mass-dependent selection efficiency determined 
from MC simulation; N{ip') = 1.06 x 10^ is the number 
oiip' events [l| and B{4j' tt+tt" J/i/') = (33.6 ±0.4)% 
is the FDG world average [l^. The upper limit is in- 
creased by a factor of 1/(1 — cr), where a is the total sys- 
tematic error (7.5%) to give a conservative result. The 
resulting B{J/ip 7^°) x B{A'^ 1^1+1.1^) ranges from 
4 X 10-'' for A° mass threshold to 2.1 x 10"^ for M{A") 
near 3.0 GeV/c^, and is shown in Fig. [Ijc). The branch- 
ing fraction upper limit is less than 10~^ for all M{A°) 
values below 0.36 GeV/c^, and is less than 10~^ for all 
masses below 2.79 GcV/c^. 
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